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Abstract. A novel determinantal representation for matrix elements of local 
spin operators between Bethe wave functions of the one-dimensional s = ^ XXZ 
model is used to calculate transition rates for dynamic spin structure factors in the 
planar regime. In a first application, high-precision numerical data are presented 
for lineshapes and band edge singularities of the in-plane (xx) two-spinon dynamic 
spin structure factor. 



1. Introduction 

The importance of the spinon quasiparticle for the understanding of the quantum 
fluctuations in integrable quantum spin chains is established on rigorous grounds 
and further supported by experiments probing the dynamical properties of quasi- 
one-dimcnsional magnetic compounds at low temperatures. Consider the familiar and 
widely studied s = h XXZ antiferromagnet, 



with J > 0, A > 0. Its ground state for even N can be configured as the physical 
vacuum for spinons and its entire spectrum can be generated via systematic creation 
of spinon pairs, for example, in the framework of the algebraic Bethe ansatz 

Theoretical and experimental evidence points to the dominance of two-spinon 
excitations in the low-temperature spin dynamics of the XXZ antiferromagnet 
|21 El EI El El • The exact two-spinon dynamic structure factor for the spin fluctuations 
perpendicular to the antiferromagnetic long-range order in the axial regime A > 1 
was recently calculated via algebraic analysis El • This method of exact analysis is 
not readily applicable in the planar regime A < 1, but exact results for the case of 
isotropic exchange (A = 1) can be inferred as a limiting case. Extending the exact 
analysis of dynamic structure factors into the planar regime has been a tantalizing 
challenge ever since. 

Remarkable advances in the determination of matrix elements from Bethe wave 
functions opened up new and promising avenues for the calculation of dynamical 
properties of the XXZ model. Of critical importance has been the recent work 
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of Kitanine, Maillet, and Terras [§] and the earlier work of Korepin and Izergin 
|1UI 1111 [T5] , which accomplished a reduction of the norms of Bethe wave functions 
and of matrix elements (form factors) for the local spin operators S", a = x,y, z, to 
determinants. We have already used their formulas to calculate explicit expressions for 
transition rates of spin fluctuation operators between eigenstates of the XXX model 
(A = 1) and have applied them to calculate lineshapes of dynamic spin structure 
factors in a magnetic field . Here we present corresponding exact expressions for 
the XXZ model in the planar regime (0 < A < 1) with an application to the in-plane 
two-spinon dynamic spin structure factor. Related projects, which focus on out-of- 
plane spin fluctuations, dimer fluctuations, XX limit (A = 0), four-spinon structure 
factors, axial regime, and T > dynamics are in progress |14| . 



2. Bethe ansatz equations 



The Bethe wave function of every eigenstate in the invariant subspace with z- 
component = N/2 — r of the total spin is specified by a set of rapidities z%, . . . , z r , 
which are solutions of the Bethe ansatz equations ^3 HB) 

N arctan ^cot ^ tanh -^j — ttJ, + arctan ( cot 7 tanh — l -— — - j , i = 1, . . . , r (2.1) 

with anisotropy parameter 

7 = arccosA (0 < 7 < tt). (2.2) 

Different solutions within this subspace are distinguished by different sets of Bethe 
quantum numbers {Ii}. It is useful to further distinguish solutions with only real 
magnon momenta ki and solutions where some or all ki are complex. The relation 
between the magnon momenta ki and the rapidities z% is 

tanh = y.i = tan ^ cot (2-3) 

The transition rate expressions to be presented below are constructed to hold for 
all Bethe ansatz solutions with real ki. The generalizations necessary to cover also 
solutions with complex ki are straightforward |17j . 

It turns out that even in the restricted set of solutions with real ki some of the 
rapidities may not be real, namely those with | tanh(zi/2)| > 1, which poses a problem 
in root-finding algorithms. To circumnavigate the problem we introduce alternate 
rapidities yi defined in (|2.3|) . They stay real for all eigenstates with real ki. The 
Bethe ansatz equations for the yi read 



N arctan (yi cot — ) = -kU + arctan ( cot 7 — — — ) , i = l,...,r. (2.4) 
V 27 j£ V ^-ViVi) 

The wave number and energy of any given solution are 



fc = 7rr -77X>> ( 2 -5) 

i=l 

g^=-y>-cosfco=-y; , sin27 , (2.6) 

J ^ v ' ^coshzi-cos7 

where Ep = NA/4 is the energy of the reference state |TT ' ' ' ?)• 



M x (g>7) = h. Vu, , 2 » M = *,+," (3-1) 
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3. Matrix elements 

In generalization to the results presented in Ref. (B, for A = 1 we use the formulas 
for matrix elements (V>o|<S'^|V'a) from Kitanine, Maillet, and Terras J2J and the norms 
from Korepin JU| of Bethe wave functions with real hi to calculate transition 

rates 

\\M 2 Hx\\ 

from the ground state of H for the operators 

^ = ^E ei9 "^' M = *,+,"■ (3-2) 

* n 

They probe the parallel (/z = z) and the perpendicular (// = +, — ) spin fluctuations 
at zero temperature. An important and largely unanticipated feature is that the 
determinantal expressions become much simpler in reciprocal space. For the parallel 
spin fluctuations our calculations yield the following results: 

N Cr({ Zi }) Q |det(H-P)| 2 

Mx{q) ~ TZ« }) ^ ({Zi}) det K({zJ) det K({*°}) ' (3 ' 3) 

where 

r r 

£r({*})=n*(*)> ^(W) =111^-^)1' ( 3 - 4 ) 

i— 1 z<j 

Hab = imhtd^z,)/^ ( n G(z ° - - d(zb) n g *(^° - *) ] . ( 3 - 5 ) 



P o6 =i2/ e («2)JjG(z i -z t ), o,6=l,...,r, (3.6) 

G(z) = sinh(z/2) cot 7 + icosh(z/2), (3.7) 
d(z) ^/tanh(,/2)cot( 7 /2)^i^ ; ^ 



\tanh(z/2) cot(7/2) 

^ J A(z a - z b ) COS7 : a 7^ 6 

i^ah — ^ Nn(z a ) — cos7^^ a K(z a — Zj) : a = b ' ^ 

k(z) = - 5 ^7 

V ' 2sinh 2 (z/2) + sin 2 ( 7 /2) V ; 

• 2 

A(z) = Sm 7 (3.11) 

sinh (2/2) + sin 7 

The corresponding results for the perpendicular spin fluctuations are 

Mf <" " " ^'W^'W' d* K^KtW}) ■ ( " 2) 
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where 



r+l r+l 



H+ r+1 = i«(z ), a = 1, . . . r + 1, 6 = 1, . . . , r, 



H * " Linh[ ( ^^)/2] I n - *> - « n ^ - -) I ■ 




H-=i/s(«2), a = l,...r, 6 = l,...,r-l. (3.14) 

In the combined limits 7 — > 0, z$ — > 0, Zi/7 — > z^, we recover term by term the 
results reported in Ref. [I3j for the XXX model (A = 1). In the XX limit (A = 0) 
considerable simplifications occur in the transition rate expressions but extreme care 
must be exercised in their applications because of singularities in the Bethe ansatz 
equations. These singularities manifest themselves, for example, in the occurrence 
of critical pairs of rapidities yi,yj in Eqs. (|2.4|) with the property lim 7 ^o UiUj = 1 
[HI USUI]. Results for the XX model are forthcoming [H] , 

More compact expressions for the results Ij3.3|l and (|3.12|) can be obtained as 
follows. Replacing (|3.9|1 by the matrix 

cos 7 K(z a - z b ) , , 

N n(z a ) r 

allows us to extract a factor N r from det K: 

det K({ Zi }) = N r C r ({ Zi \) det K({zJ). (3.16) 
Using the Bethe ansatz equations in the algebraic form 

again in H a b, P a b, and H^ fc , further simplifies these matrices. The consolidated 
expressions then read 

= N JC r ({z°}) |det(r-fl)| 2 
x[q> 4 K r {{ Zi }) detK(K})detK({*})' ^ j 

,.+ , , _ f lCr±i{{z l }) \ ±1 Idetr^l 2 

" I XV(K}) J det K({*}) det K({z«}) ' ^ iJJ 

where 

rab ^ Fn[z - Zb) n G(^y + G^^y n ^-^y J • (3 - 20) 

r + j_ jp / o\ I G ( z r+i - Zb) tt G{z - zg) G*(z r+ i - zg) -j^ G*(zj - zg) 
ab ^ [ G(z a - zg) Ij G(z° - zg) + G*(z a - zg) 11 G*(z° - zg) j ' 

r+ r+1 =l, 0= l,...r + l, b= l,...,r, (3.21) 
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G(z° - z b ) 1 = 1 G(z, - z b ) ' G*(*° - z 6 ) 1 = 1 G*( 2j - z 6 ) 
r-=l, a=l,...r, 6=l,...,r-l, (3.22) 
P , /\ • sin 2 ( 7 /2)+ S inh 2 (z/2) 

Fw(z ' z) "iVsin 7 sinh((,-,')/2)' ( } 

/ 1 ••• 1 



\ i ... i 

In the isotropic limit 7 — ► we can use the results with 

and the other ingredients as defined in |13| . 

The corresponding results for M?(q) and M^(q) expressed in terms of the 
rapidities yi are equally compact |14j . For applications involving excitations with 
real magnon momenta ki, we solve the Bethe ansatz equations l|2.4l) and calculate 
transition rates from ({3.18(1 and H3.19[l . 

4. Two-spinon dynamic structure factor 

To demonstrate the practicality and usefulness of the newly derived transition rate 
expressions, we present high-precision data for the two-spinon part of the T = 
dynamic structure factor 

S^ + (q,u>) = 2irY,M;(q)6(u J -w x ), (4.1) 

A 

which captures the in-plane spin fluctuations. The state \ipo) in ijOJ is the XXZ 
ground state (spinon vacuum). It has quantum number = and is characterized 
by the solutions of Eqs. (|2.4|l for the set of r — N/2 Bethe quantum numbers 
n ' = — (N + 2)/4 + i, i = 1, . . . , N/2. Excited states containing two spinons with 
parallel spins are known to account for a major contribution to the in-plane spin 
fluctuations 

The number of two-spinon states with 5f, = 1 is ^N(N + 2). Their r = N/2 - 1 
Bethe quantum numbers, which are integers for odd r and half-integers for even r, 
comprise all configurations 

N N 
- -j < h < h < ■■■ < It < -j. (4.2) 

These states are described by real solutions of Eqs. I|2.4|l . For N — > 00 they form a 
continuum in (q, w)-space with boundaries |15ll3"| 

, . 7rJsin7 ... , . 7rJsin7 . q 

e M<?) = — z sing , eu(q) = sin- . (4.3) 

z 7 7 z 

The scaled density of two-spinon states constructed from D^> (g, uj\) = 2tt/[N(uj\+i — 
u>\)] for the energy-sorted sequence of two-spinon excitations at fixed q and finite N 
turns into the function 

D^{q^) = [el{q)-^]-^ (4.4) 
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for clW) < lo < €jj{q) m the limit N — > oo. For the two-spinon part of the dynamic 
spin structure factor we use the product representation 

Si a i(«,w) = M^l(q,uj)D^( q ,cj) (4.5) 
as discussed in a previous application of a similar nature |21| . The first factor 
represents the scaled transition rates M_+(q, lo) = N\ (G\S~ |A) | 2 between the ground 
state and the two-spinon states (|4.2() . 

( 2 ) 

Here we focus on the lineshape and on the singularity structure of S_i_(q,ui) at 
q = 7r, where the spectral threshold is at lo = 0. In Fig.^ we show N — 512 data for 

(2) 

the transition rate function Mlj_(7r,u;) at A = 0.1,0.3, . . . ,0.9. The data indicate an 
infrared divergence and exhibit a monotonically decreasing w-dependence toward zero 
intensity at the upper band edge. The inset to Fig. ^ zooms into the behavior near 
eu{^)- The data points approach zero linearly with a slope that becomes smaller with 
decreasing A. A linear behavior was already established via algebraic analysis for the 
case A = 1 [7]. 




CO/J 



Figure 1. Scaled transition rates between the XXZ ground state and the two- 
spinon states at q = it for N = 512 and A = 0.1, 0.3, . . . , 0.9. The inset shows the 
same data again for ui near ejj{n) on a different scale. 

The spectral- weight distribution of the two-spinon contribution to S y(n,Lu) is 

then inferred via product representation (|4.5|l from the transition rate data and the 
two-spinon density of states 14.4(1 . The results are shown in Fig. The divergent 
density of states ((4.4(1 at u = ec/M converts the linear cusp of the function M^(ir, lo) 

(2) 

into a square-root cusp in the function S_+(tt,lo) as is best visible in the inset. The 
A-dependent exponent of the power-law singularity is exactly known (221 12.'{| : 

S-+(tt,lo) - lo^ 1 ^^. (4.6) 

This non-universal critical singularity is accurately reflected by the N = 512 data as 
is demonstrated by a data fit at low frequencies for each lineshape. 

In summary, the availability of determinantal transition rate formulas opens up 
a whole new area of applications of the Bethe ansatz with enormous potential for 
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Figure 2. Lineshape at q = it of the two-spinon contribution to S y(q,w). The 

dots represent data for N = 512 at A = 0.1, 0.3, . . . , 0.9 (left to right). The lines 
are two-parameter fits aw -1 1/ v + b (over the interval shown) with the exactly 
known exponent from Eq. 14.61 for N = oo. The inset shows the N = 512 data 
again for lo near e[/(-7r) on a different scale. 

important new results, among them results of relevance for experiments (neutron 
scattering, NMR) on quasi-one-dimensional magnetic insulators. 
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